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Abstract 

For two-dimensional lattice equations one definition of integrability is that the 
model can be naturally and consistently extended to three dimensions, i.e., that it 
is "consistent around a cube" (CAC). As a consequence of CAC one can construct a 
Lax pair for the model. Recently Adler, Bobenko and Suris conducted a search based 
on this principle and certain additional assumptions. One of those assumptions was 
the "tetrahedron property", which is satisfied by most known equations. We present 
here one lattice equation that satisfies the consistency condition but does not have 
the tetrahedron property. Its Lax pair is also presented and some basic properties 
discussed. 

1 Introduction 

Within the field of integrable dynamics there have recently been interesting developments 
in the study of integrable difference equations (an overview of the topic can be obtained 
from the proceedings mentioned in pQ). These include both the discussion on what is a 
proper definition of integrability, or whether the different suggested definitions actually 
agree or not. Numerous integrable (in some sense) difference equations have been proposed 
by discretizing known ODE's and PDE's in a way that retains some good properties. Once 
a definition of integrability has been proposed one can also try to search for all equations 
having the chosen property. 

In this paper we consider integrable difference equations denned on a 2-dimensional 
lattice. We assume that the lattice is rectangular and infinite, and concentrate on some 
square in it (see Figure GJ (We normally use subscripts in square brackets to indicate 
shifts in the indices, but for composite expression we sometimes also use ~ for shift in the 
1-direction and ~ for a shift in the 2-direction.) In general the map is given in terms of a 
multi-linear equation relating the four corner values in Figure GJ 

K xx[i\X[2]X[ 12 ] + L\ xx^xp] + L 2 xx[ 1 ]X[ 12 ] + L 3 xx[ 2 ]X[ 12 ] + L4 x^x^x^] 

+Pl XX[l] + P 2 X[i]X[ 2 ] + P3 X[ 2 ]X[ 12 ] + P 4 X[ 12 ]X + P 5 XX[ 2 ] + P 6 X[!]X[ 12 ] 

+QlX + Q 2 X[i] + Q 3 X [2] + Q 4 X[i2] + U = Qi 2 (z,qi],Z[2],£[i2];Pi,P2) = 0. (1) 
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Figure 1: The lattice map is defined on a elementary square of the lattice 

Here the coefficients K, L u , P u , Q u , U may depend on the two spectral parameters Pi,P2- 
If any 3 of the corner values are given then the fourth one can be obtained as a rational 
expression of the other three. One can therefore propagate any staircase-like initial value 
line to cover the whole plane [2 . 

How should integrability be defined for such maps? In [Sj the following definition 
of integrability ("Consistency Around the Cube", CAC) was proposed (see Figure |2J): 
Adjoin a third direction (therefore assuming x = x n ^ m ^) and use the same map (but with 
different spectral parameters) also in planes corresponding to indices 1,3 and 2,3. That is, 
the map given in contains shifts and parameters associated with directions 1,2 now the 
same should be done with directions 3,1 and 2,3, furthermore, on on the parallel shifted 
planes we use identical maps. We assume that the values x, x^j, xp], £[3] at black circles 
in Figure 121 are given, then the values at open circles are uniquely determined using the 
relevant map, but the value at #[123] can be computed in 3 different ways, and they must 
give the same result. In other words: 
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£[12] 
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£[i], 


X[ 2 ], 


£[12] 


Pi, Pa) 


= 0, 


solve 


£[23] 
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= 0, 
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Figure 2: Given the values at the black circles, one should get a unique value for £[123], 
even though there are three possible way to compute it. 
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then 2?ri23] computed from 



Q 12 (X [ 3 ] , X [31] , X [23] ,X[i23]',Pl,P2) 
Q23{X[l] , £[12] , £[31] , £[123] ] P2,Pz) 
Q31 0&[2] , 3J[23] , &[12] , aJ[i23] ; P3, Pi) 







0. 
0. 



or 



or 



should be the same. The functions Qij could be in principle be different, but are usually 
assumed to be identical. 

The above idea (and diagram) resembles, e.g., the 3D Bianchi diagram that is obtained 
from consistency of Moutard transformations (see j3J), but it is here used in a different 
context: the diagram introduces a condition for lattice maps defined on squares rather 
than for Moutard transformations defined on edges. Furthermore, this CAC principle is a 
constructive definition of integrability in the sense that it leads algorithmically to a Lax 
pair. 

The following maps are well known examples that have CAC property PJ Ej: 

1. Lattice KdV: (pi - p 2 + £p] - £[i])(pi + Pi + x - X[i 2 }) = p\—p\, 

2. Lattice MKdV: pi(xx[ 2 ] - £[i]£[i2]) = M^pL] - X[2]X[i2]), 

3. Lattice SKdV: {x - aj[ 2 ])(2C[i] - ^[i2])p| = (x - X[i])(x[ 2 ] - X\i2])p\- 

2 Searching for integrable lattices 

Now that the definition of integrability has been given one may ask for a listing of all 
integrable models. The complete classification of CAC maps is in fact a formidable open 
problem. The set of equations that needs to be solved can be derived from above: compar- 
ing the different forms for X]vm and collecting the various coefficients of x and its shifts, 
we get 2 x 375 functional equations for the coefficient functions (in practice we have often 
used 3 x 375 equations for a symmetric approach). The equations are polynomial in the 
coefficient functions appearing in (JTJ), but the functions depend on different pairs of the 
three spectral parameters A,, i = 1,2,3 and all three Aj's appear in each equation. The 
equation list itself takes 35MB to store. 

In |Hj the consistency equations were solved in the case where the map Q is the same 
on all planes, and under two additional assumptions: 

• symmetry: 



• "tetrahedron property" : £[123] does not depend on x. 

With these assumptions the authors were able to get a full classification resulting with 9 
models. 

The tetrahedron assumption is indeed satisfied by most of the well known models, but 
one can nevertheless ask whether it is a fundamental or essential property and whether 



Q(x, x m , x[2], X[ 12 ])Pi,p2) = eQ(x, X[ 2 ], X[x], X[i 2 ];p 2 ,p 1 ) 
= o-Q(x [1] ,x,x [12] ,x [2] ;p 1 ,p 2 ), e,a = ±1, 



(2) 
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there are any nontrivial models that do not satisfy it. It should be immediately observed 
that the more or less trivial models 



X[i]X[ 2 ] - ££[i2] = 0, and y - y m - y [2 ] + y[ 12 ] = 0, (3) 

(x = e y ) do have CAC property but not the tetrahedron property, in fact one quickly 
finds that 

£[123] = £[i]£[2]£[3]£~ 2 , and y {12 3] = y\i] + yp] + V[3] ~ 2y, (4) 

respectively. One may now ask whether these models have nontrivial extensions. 

The symmetry assumptions (J2J) lead to 2 possibilities w.r.t. a, the one with a = — 1 is: 

Q = ai(pi,p 2 )(££[i]£[2] - ££[2]£[i2] - ££[i]£[i2] + £[i]£[2]£[i2]) 

+a 2 (pi,p 2 )(xx [12] - £[ipp]) + a 3 (p 1 ,p 2 ){x - £[ij - x [2] + x [12] ). (5) 

In |5] it was observed that in this symmetry class there are no integrable cases with 
tetrahedron property, which is true. The ansatz (jSJ) forms, fortunately, a rather simple 
class and a direct computation shows that it contains one new integrable map 

££[i]£[2] - ££[1]£[12]-££[12]£[2] + £[i]£[2]£[i2] 

+ (£[i]£[2] -££[i2])(Pi +P2) -P\P2{x ~ - %[2] +£[12]) = 0, (6) 

with two parameters. Since it can also be written as 

£ £[i] £[2] -pq £[12] 

-£ £[ij £[i2] + (p + q) £[i] £[ 2 ] +pq £[2] _ 
-£ £[2] £[12] ~(p + q) x £[12] +pq £[i] 

+£[1] £[ 2 ] £[12] ~pq X 

we call it "the bow-tie model" . Model © is, however, still Mobius-equivalent to the first 
model in © by x — > — (p\X +p 2 ) /{x + 1) [although this transformation does not trivialize 
the maps on the other planes, since it depends explicitly on pi-,p 2 }. Therefore we decided 
to search for possible extensions. This was done perturbatively, with completely general 
1st and 2nd order extensions to the bow-tie map (this work was done with REDUCE 
3.70, and 1GB memory). In this approach it is not necessary to specify the details 
of the parameter dependence - indeed the result was surprising in this respect. From 
the outcome of this exercise we noticed certain weaker symmetry properties among the 
coefficient functions and when the (nonperturbative) computations were done with these 
properties (and translational invariance) the following solution was found: 

Q(x, £[i] , £[ 2 ] , £[i2]; ei, oi;e 2 ,o 2 ) = 

££[l]£[ 2 ](Oi - 2 ) - ££[l]£[l2](ei - 2 )-££[ 2 ]£[12] («1 ~ e 2 ) + £[1]£[2]£[12] ~ e 2 ) 

+ (x[i]X [2] - ££[i2])(eiOi - e 2 o 2 ) 
+{xx[ 2 ]Oi - £[i]£[i2]ei)(e 2 - o 2 ) + (£[2]£[i2]e 2 - ££[i]0 2 )(ei - oi) 
-£[i2]eie 2 (oi - o 2 ) + £[ 2 ]e 2 0i(ei - o 2 ) + £[i]ei0 2 (oi - e 2 ) - £Oi0 2 (ei - e 2 ) 

= 0. (7) 

This new model is the main result of this paper. It is interesting to note that in (J7J) there 
are two parameters Oj,ej in each direction, and that if e, = Oj then 0\ — o 2 factors out 
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leaving (jfij). The parameters with different indices must all be different, otherwise the 
map factorizes. This model obeys the CAC property, in fact one finds 

num(a;[i23]) = -2?(z[i] + 0i)(x[ 2 ] + o 2 )(x[ 3 ] + o 3 )(oi - o 2 )(o 2 - o 3 )(o 3 - oi) 
+ (x[i] + oi)(x [2] + o 2 )(x [3] + o 3 ) 

x [(eie 2 + e 3 o 3 )o 3 (oi - o 2 ) + (e 2 e 3 + eio x )oi(o 2 - 03) + (e 3 ei + e 2 o 2 )o 2 (o 3 - 01)] 
+ (a; + e 3 )(a;[i] + oi)(x[ 2 ] + o 2 )o 3 (oi - o 2 )(e 2 - o 3 )(o 3 - e x ) 
+ (x + ei)(x [2] + o 2 )(x [3 ] + o 3 )oi(o 2 - o 3 )(e 3 - oi)(oi - e 2 ) 
+ (x + e 2 )(x [3] + o 3 )(x[i] + oi)o 2 (o 3 - oi)(ei - o 2 )(o 2 - e 3 ), 
den(x[i 23] ) = (sqi] + oi)(x [2] + o 2 )(x[ 3 ] + o 3 ) 

x [(eie 2 + e 3 o 3 )(o 2 - 01) + (e 2 e 3 + eiOi)(o 3 - o 2 ) + (e 3 ei + e 2 o 2 )(oi - o 3 )] 
+ (x + e 3 )(x[ij + oi)(x[ 2 ] + o 2 )(oi - o 2 )(ei - o 3 )(e 2 - 03) 
+ (a; + ei)(x [2] + o 2 )(x [3] + o 3 )(o 2 - o 3 )(e 2 - oi)(e 3 - 01) 
+ (x + e 2 )(x [3] + o 3 )(x [1] + Oi)(o 3 - o x )(e 3 - o 2 )(ei - o 2 ), 

which are symmetric under permutations of 1, 2, 3, and the explicit x-dependence demon- 
strates violation of the tetrahedron property of [5]. 



3 Symmetries 

The model (J2J) is invariant under the Mobius 

aX + /3 
jX + 5' 

It has the symmetries 

Q(x,X[i],X[ 2 ],X[i 2] ; ei, o 1( e 2 ,o 2 ) = 
Q(x, , xp], X[ 12 ] ; ei, 01, e 2 , o 2 ) = 



transformation 

X, X2, X^', X^'fc 

P = ei or Oj, 

-Q(x, X[ 2 ], X[i] , X[i 2] ; e 2 , o 2 , e 1( o x ), (8) 
-Q(x[i],x, X[i 2] , X[2]; 01, ei, e 2 , o 2 ). (9) 



These are similar to (0) with e = er = — 1, but note that in Q there is an exchange 
between the two parameters related to direction 1. Indeed, it seems that the extension 
of the parameter space from 1 to 2 dimensional (in each direction of the lattice) allows 
the introduction of a nontrivial exchange symmetry, resulting with the nontrivial model. 
(Recall that if = o, the model simplifies to 

The map (JJJ) can also be written in the following rational form 

x + e 2 X[i2] + 02 _ X[!] + e 2 x [2 ] + o 2 
x + e\ X[i 2 ] + oi X[i] + ox X[ 2 ] + e\ 

(The Backlund equation presented in Eq. (5.11) of [Zj is similar, but the distribution of 
the terms is in fact different.) A form equivalent to (jl(lj) is 

02 + X[ 12 ] e 2 + x _ oi + X[ i2 ] e 1 + x 
o 2 + x [2 ] e 2 + X[ij ox + X[i] e 1 + x [2 ] ' 



5 



and comparing these one observes duality under 



o 2 <-> x, e 2 <-> x [12 ], ox <-> ei <-> xpj. (11) 
We also note that (jZJ) has the simple solution 

Xn,m, P = c + (ei - oi)n + (e 2 - o 2 )m + (e 3 - o 3 )p. 

4 Lax pair 

An important property of the CAC definition of integrability is that it provides a method 
of constructing a Lax pair. The recipe was provided by Nijhoff in [H]: We consider the 
three maps on planes 12, 23, 31. The third direction is taken as auxiliary (spectral) and 
the system is linearized in the corresponding variable X[ 3 ] and its shifts. We also use 
notation o 3 = /i, e% — A. 

Solving for x [31] from Q(x, x [3] , arpj , £[31] ; A, fi, e t , 01) = yields 

X[ 3 ][aa?[i](oi - n) +xo 1 (\ -n) + X[i](e x oi - A//) + Aoi(ei - y)\ 
+xx[i]n(—e 1 + 01) + xoi/i(-ei + A) + X[ijei/i(-A + 01) 

[31] x [3 ] (-ei + A) + x(-X + 01) + A(-ei + 01)] 

(ei — //) + x(eiOi — A/x) + X[ijei(A — /i) + eiA(oi — ji) 

and this is then linearized by introducing /, g by 

/ /[2] /[I] 

Z[3] = -, X[23] = — , *[31] = — , (12) 
# #[2] P[l] 

resulting in 

/[1] =«i / [^[i](/i - 01) + xoi(/i - A) + X[ij(A/i - eioi) + Aoi(/i - e x )] 
+Ki g [xx[i]fi(ei - 01) + xoin(ei - A) + X[ijei/i(A - 01)] , 
=«i / [z[i](ei - A) + x(A - o x ) + A(ei - 01)] 

[xx - ei) + x(A// - eioi) + cque^// - A) + eiA(/i - 01)] . 

Here the overall separation factor /«, may contain x, xuy We write this, and the corre- 
sponding equation obtained from Q(x, X[ 2 ], a; [3], £[23] j e 2, o 2 , A, /x) = in matrix form, 

with 



U = Ki(x K + Oi)(x + A) ^ ^ J + K,(a; H + A) (a; + e { ) ^ ^ ^ j . (14) 

Now from the consistency condition (0[i])r 2 i = (0[2])m we get the matrix equation 

L\L 2 = L 2 L U (15) 
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which is satisfied modulo the map ( | 7 j) . provided that the separation constants Kj satisfy 



Ki«2 _ (gjgj + giXg + e 2 ){x [2 ] + A) 
KiK 2 + ei)(x[i] + e 2 )(x[i] + A)' 

A simple solution to this is 

1 



(x + 6i)(x + A) 



(16) 



(17) 



leading to 

^ l x + 6i \-l —e%) x + A V 1 /i 
A similarity transformation with S = ( -i ° ) simplifies this further to 



X[i\ + Oi ( j_M X[i] + A f—Oi —fio, 



( XM+Oi XM+Oi X[i\+X\ 

I [V ~ Zi) ; m : m — 

X + Ci X + 6i X + A 

V (/i -° l) TTX / 



(19) 



5 Summary 



The "consistency around the cube" definition of integrability is very transparent since 
it directly leads to the Lax pair. It also defines a clear search problem whose complete 
solution is, unfortunately, still beyond our computational abilities. The integrable model 
(jZJ) or fllOl) presented in this paper is the first result in our search project. This model has 
several interesting features: 

• It does not satisfy the tetrahedron property. 

• It has two parameters in each direction. 

• It is dual under interchange of variables and parameters. 

Further properties of this model will be discussed elsewhere. 

In this general class of lattice models (multi-linear, consistent around a cube) there 
are probably still many other models to be found. 
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